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I. INTRODUCTION 

Cosmic strings are now a widely recognized part of cosmo- 
logical theory. Cosmic stings appear naturally in a multitude 
of inflationary models as topological defects from the early 
Universe (e.g. [UE], for more studies see references in [3]). 
Similar objects commonly referred to as cosmic superstrings 
can also be produced in fundametal string and M-theories 
[4] [5] . In the present study we will not differentiate between 
the two classes, because their observational signatures consid- 
ered in this paper are the same. 

The key parameter of a cosmic string is its tension /i, which 
is assumed to be related to the effective energy scale of the 
string-producing theory A by [H] 
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Earliest theories of string formation placed time of their gen- 
eration to the Grand Unification Theory epoch and there- 
fore their tension seemed to be of order 10 -6 . Initially, 
possible tensions of string were constrained from both sides: 
10 -11 < G/i/c 2 < 10~ 6 , but eventually the lower bound was 
removed and strings with arbitrarily low tension are theoret- 
ically allowed now (e.g. [7]). Cosmic strings with low tension 
can solve some astrophysical problems: e.g., recently, cosmic 
strings with tensions about Gfi/c 2 ~ 10~ 12 were proposed as 
prominent source of high-energy cosmic rays [8] . 

Simulations [T] [9] suggest the energy fraction in strings has 
a scaling behavior: their density Q s (in units of critical den- 
sity po = 3Hq /8iyG) does not depend on cosmological time. 
Recent works suggest that strings contribute a subdominant 
fraction to the energy balance of the Universe [TO] . This ques- 
tion has not been completely resolved yet. We will treat f2 3 
essentially as a free parameter and will try to constrain it 
observationally. 
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There have been a number of attempts to limit the den- 
sity of strings from various observational perspectives. Heavy 
enough strings, if present, would make distinctive imprints 
on the cosmic microwave background and shapes of dis- 
tant lensed galaxies; these techniques indicate absence of 
strings with G/i/c 2 > 10~ 7 in the Universe [UJ HI fT3j . 
More stringent, though more model-dependent constraints 
come from pulsar timing (PT): cosmic strings emit grav- 
itational waves and the corresponding background can be 
detected by usual methods of PT[14|. This method rules 
out any significant presence of cosmic strings with tensions 
fi at Gfj,/c 2 < 10" 9 - 10" 8 [IS]. Ultimately, planned mis- 
sion LISA are expected to test the presence of gravitational 
wave background from lighter strings with tensions down to 
Gfi/c 2 ~ lO" 14 - 10- 16 . 

String network has a complicated structure with a combi- 
nation of long straight segments and a population of string 
loops of various lengths L that were formed in the intercon- 
nections between straight strings. The loops oscillate rela- 
tivistically (f3 ~ 0(1)) with an amplitude of order L and pe- 
riod T = L/2c, emitting gravitational waves and eventually 
decaying. Only sufficiently long loops survive by the present 
time. 

Recent simulations 3 show that the surviving large-scale 
(L ~ 1 pc and above) loops of light cosmic strings experience 
considerable clustering, which closely follows that of the dark 
matter, albeit with a somewhat lower amplitude. In central 
parts of large galaxies, such as the Milky Way, the loops' 
density can be enhanced by up to 10 relative to its average 
cosmological value. This density enhancement significantly 
boosts the detection rates of experiments sensitive to the local 
population of lenses [31 [16] . 

In this paper, we investigate how the clustering affects the 
prospects of detection of local cosmic strings via two obser- 
vational signatures - lensing on the string as seen in the pho- 
tometry of background objects and Kaiser-Stebbins effect af- 
fecting the timing of millisecond pulsars. 

The paper is organized as follows. In the next Section |TT| 
we discuss the expected rate of the effects and lay out a simple 
formalism with which one can interpret the non-detection of 
the effect in a given experiment. Then in Sections |III| and |IV| 
we consider lensing and Kaiser-Stebbins effect in a greater de- 
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tail and calculate the constraints on the string loop population 
from existing observations. The final Section discusses the ob- 
tained results and offers suggestions for further research. 



II. EVENT RATE AND CONSTRAINTS FROM 
NON-DETECTION 

The first effect we consider is lensing of a background object 
by the string, in which two identical positive-parity images of 
the object appear on the sky when it enters a narrow strip 
along the string (e.g., |f 71 If 8|). The width 8 of the strip 
depends on the tension fi of the string and the local inclination 
8 of the string to the line of sight: 
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where D os and D\ s are the distances, respectively, from the 
observer to the source and from the string to the source (along 
the line of sight). For the presently allowed tensions of the 
string, the two images cannot be resolved and only an increase 
in the total brightness of both images is seen as the source 
traverses the strip. 

The second observational signature of our interest is the 
Kaiser- Stebbins effect, which stems from the induced Doppler 
shift in the conic space-time around a moving string and re- 
sults in a discontinuity A/ of the observed frequency / of any 
radiation from the source observed on the either side of the 
string [1SJ 120] • It is a l so proportional to the string tension 
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where 7 = (1 — /3 2 )~ 1,/2 and f3± are, respectively, the Lorentz- 
factor of the string and the orthogonal (to the string) compo- 
nent of its transverse (to the line of sight) velocity w.r.t. the 
source (in units of c). By itself, this frequency jump is again 
too low to be observed directly (e.g., spectroscopically) but 
it can manifest itself in the integrated form of pulsar timing 
residuals. 

Both the lensing and Kaiser-Stebbins effects depend on the 
crossing of an observed source by a cosmic string. Thus, the 
probability that a given source if affected by either effect dif- 
fers from the standard value of po = f — exp(— ro) where to 
is the optical depth given by the fraction of the sky covered 
by cosmic strings. However, to is still useful for an estimate 
of the (un) importance of strips overlap in lensing. 

For cosmic strings to is given by summing the contribu- 
tions of all infinitesimal slices of thickness AD within a solid 
angle d 2 <9 along the line of sight D S (0,D OS ). Each contri- 
bution dTo is the fraction of area covered by all strips within 
the volume D 2 d 2 0dD relative to the total area of this slice 
D 2 d 2 0. The area covered by all strips within the slice is given 
by their total length p s j jj,D 2 dD times the linear width DS of 
strips: 

dT = (DVer 1 ^D 2 d 2 6dD xD^ —DSdD. (3) 

Given that 8 ~ p (Eq[T]), r is actually independent of fi: 
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Equation Q is the same as the optical depth due to a pop- 
ulation of point lenses with mass density 2p(|sin#|) = np/2. 
These values have been estimated in a number of works both 
for the local case of MACHOs in our Galaxy and for a hypo- 
thetical cosmological population of compact lenses. As long 
as the density in strings does not exceed that of the dark mat- 
ter, we can use the theoretical upper limits from those works 
to constrain to. The estimates (to < fO -6 — fO -5 for the 
Galaxy and to < fCF 2 — fCF 1 in the cosmological case with 
source at redshift f '21, 22, 23, 24 ) are significantly below 
unity. This allows us to assume that the strips do not signifi- 
cantly overlap in projection and at every given moment every 
source is affected by just one string at most. 

However, the optical depth defined by |3| is not a measure 
of the expected rate of events in an observational search for 
cosmic strings. The latter is given by how often a given source 
is crossed by a cosmic string and given that by their very 
nature strings are extended objects, as they move across the 
sky they quickly sweep areas much greater than those swept 
by Einstein circles around point lenses with the same optical 
depth and velocities. 

Every source swept by a string will split into two parts sep- 
arated by |T]) and their combined brightness can jump by up 
to a factor of 2 while the observed frequency of the source will 
experience a discontinuity of (J2J . Whether (and how) these 
effects can be detected is a separate question depending on 
the details of a particular observational survey and discussed 
in the following. However, when estimating the optical depth 
similarly to (J3j> , instead of SD i width one should use the 
width of the strip swept by the cosmic string over the time 
span T of the observational survey; this is given by c{j3±)T, 
where {/3±) is some typical value of fix among strings con- 
tributing to the optical depth. The correct equation for the 
latter is then 1 
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which can be made very large by lowering p, due to the very 
poor existing lower limits on the string tension. The prob- 
ability p that a change - in either brightness of frequency 
of the observed source - occurs in a given source during the 
observations is p — 1 — exp(— t). 

The model of the experiment we consider is an observa- 
tional survey searching for either a characteristic transient 
brightening of the source or a jump in its frequency. The time 
span of observations is T and the number of sources monitored 
is N. The cosmic strings population is characterized by its 
average cosmological density fl s in critical units and tension 
fj, (typical length of loops L cannot be constrained with ob- 
servations we consider). We assume that the string loops can 
cluster as suggested in [3] with a local enhancement r\ over 
their average cosmological density f2 s po so that their local 



where we used the average (|sin(9j) = 7r/4 and added a sub- 
script 'ol' to D to reflect its role in traditional lensing. 



1 Strictly speaking, one cannot simply add up infinitesimal contri- 
butions because nearby layers dD are not strictly independent — 
that is, two neighboring layers most likely either will both have 
a piece of a string in them or neither will have one. However, 
assuming that the length of the string is small compared to the 
distance to the source, one can choose a 'physically infinitesimal' 
depth of the layer dD that will make layers effectively indepen- 
dent. 
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where Hq is the Hubble constant; we use rj from [3j. 

Let us now use Bayesian inference to see how a non- 
detection of the effect in the survey constrains the parameters 
of the strings (fi, p). Given these parameters, we can calculate 
the probability that the effect happens in i-th source: 
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integrating the enhancement r\ along the line of sight. 

The probability of actually detecting the effect needs to 
account for the limited efficiency of any survey. Even if the 
change in brightness or frequency does occur, we might miss 
it because the observations are not continuous, the effect is 
too weak to be observed or can hide in the intrinsic variations 
of the source. This is accomplished by multiplying pi by ef- 
ficiency factors yet that measure the probability of detection 
given that the effect does happen. 

Assuming that the sources are affected by strings indepen- 
dent of each other, the probability Q of non-detection in the 
entire survey is then the product of the probabilities of non- 
detection in every source 

N 

Q = V (No detection|n s , \i) = Y\ (1 — XiPi) , 

8=1 

which, in the Bayesian sense, evaluates to likelihood 



Q(n a ,fi) = H\ i 
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The posterior distribution density of (Q, s ,p) set by the non- 
detection is then given by the Bayes theorem: 



p (fi s , p) = Q(Q S , p)p(il 3 , p) 
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where integration in the denominator extends over the pa- 
rameter space measured by the prior p(Q 3 ,p). 

Alternatively, one can interpret p'(Q s ,p) as a scaled prob- 
ability density of Q 3 for a fixed p so that the probability of 
string density being less than Q as a function of p is 



P(SI.,H) = j P '(n' 3 ,p)dQ' s J 



p'(Q' s ,p)dfl' 3 



(8) 



Efficiencies depend not only on (Q 3 ,p) but also on the 
properties of the survey, including those of individual sources. 
They will be calculated below for both lensing and Kaiser- 
Stebbins effects. 

A choice of Xj G {0, 1} is particularly convenient. Assum- 
ing that rj is constant within the probed volume of space, the 
likelihood function reduces to a simple form: 



Q(fl s , p) = exp 
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where the sum includes only those targets, for which h% = 1. 

For this likelihood equation (JsJ> simplifies significantly once 
a flat prior in il s is chosen. In this case Q(£l s ,p) = 
exp [—Q s /Sl(p)], where the sensitivity of the experiment 
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This value then constrains the string density Q 3 at a given 
confidence level P according to 



n 3 (p,P) = -Slfjj,) ln(l-P) 

Gp 87rclnl/(l - P) 

These equations will be used below. 
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III. LENSING BY COSMIC STRINGS 
A. Phenomenology of string lensing 

Cosmic strings produce a distinct pattern of lensing 
(e.g., [171 118]). For every point-like source entering a nar- 
row strip of width 5 along the string, a second positive-parity 
image appears in a 'duplicate' strip on the other side of the 
string, the two images separated by 8. As the source contin- 
ues its way towards the string, the second image moves away 
from it until the first image disappears on contact with the 
string itself (see upper plots in Figure [T]). 

If the two images cannot be resolved, which is the case we 
consider in this paper, neither an astrometric shift in the po- 
sition of the source on the sky can be detected, and what the 
observer sees is a temporary increase in the total brightness of 
the source. For a small enough source that fits into the strip 
completely, the brightness will increase twofold. For larger 
sources, the increase will be given by the total flux in the 
part of the source that fits within the strip. 

In particular, for a disc of uniform brightness, simple ge- 
ometry gives the following fractional increase in the area and 
flux: 
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where r is the radius of the source and x is the distance from 
its center to the mid-line of the strip. 

For small sources, r ^ 5/2 the amplitude of the effect 
/max = 1 while for r > 8/2 
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For very large sources r > iS the light curve and amplitude 
formulae reduce to 
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respectively. Figure [T] shows the morphology of lensing and 
the respective light curves for uniform discs of different radii. 

Light curves for more complex sources would be similar and 
for sources smaller than the strip width the curve will be the 
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FIG. 1: The morphology of images and expected light curves 
for cosmic strings lensing a uniform brightness disc of various 
sizes. From upper-left to bottom-right r/S is 0.2, 0.5, 2, 5. The 
position of the string and the strip edges are shown by solid 
and dotted straight lines, respectively. 



same as in the top left panel of Fig. [T] except for the precise 
shape of the 'wings' leading up to the 'plateau'. Most impor- 
tantly, the amplitude - i.e., height of the plateau - will be the 
same, / max = 1, and this is a distinct photometrical signature 
for string lensing of small sources; we are not aware of any 
other phenomena that produce an exactly twofold increase in 
brightness naturally. 

The time scale t of such events is given by the angular 
velocity {f3±)c/D \ of the string with respect to the source in 
the plane of the sky and is given by 
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where 2r = 2R/D os now stands for the overall size of the 
source in the sky. 

In summary, one can see that two different regimes of lens- 
ing are possible, depending on the size of the source. For large 
sources r 3> 5, the amplitude and time scale of the effect de- 
pend on the size of the source 
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For small sources r < i, the amplitude saturates while the 
time scale depends mostly on the string tension: 

_ 2 G\x A>iAs 
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For a rough estimate as to which regime applies based on the 
observed properties of a thermally radiating uniform source, 
one can utilize the Stefan-Boltzmann law to see that the an- 
gular size depen ds on the observed flux and temperature only: 
r = T~ 2 \J F I o . More accurate estimates need to account for 
the interstellar absorption and for the actual spectral energy 
distribution and its overlap with the detector's bandwidth. In 
astronomical terms, one can obtain 
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where m is the apparent instrumental magnitude of the ob- 
ject, T e ff is its effective temperature, A is the interstellar ab- 
sorption and BC is the bolometric correction. Both T e g and 
BC can be estimated from the apparent color of the object 
while A is well mapped in the sky. 



B. Optical photometry: CoRoT and Kepler 

For hydrogen-burning, Main Sequence stars which make 
up the majority of the Galaxy's stellar population, the linear 
size R varies within roughly (10 10 — 10 12 ) cm and at distances 
D of order 1 pc to 30 kpc this corresponds to angular size 
r € (10 -13 , 10 -7 ); including white dwarf and red giant stars 
widens this interval by about an order of magnitude on both 



ends. Therefore, for strings with G/x/c 
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small-source and large-source regimes may be relevant in the 
case of the Galaxy, though numbers tend to point towards 
the latter case. The amplitude of the effect could be any- 
where from 10 -7 to 1 while the time scales of cosmic string 
lensing events within the Galaxy are rather short. Even in the 
extreme case, where a red giant of size R ~ 10 13 cm seen at 
the far end of the Galaxy (D ~ 30 kpc) is crossed by a string 
of tension G/i/c 2 ~ 10~ 10 located half-way to the source, ap- 
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plication of ( 13 1 yields t of order an hour. For a more typical 
case of a main-sequence star at a few hundred pc from the 
Sun this reduces to tens of seconds or below. 

There is some overlap of these estimates with the target 
sensitivity of exoplanet-hunting missions such as CoRo T and 
Kepler. Their target stars are relatively close to the Sun 
(~ 100 pc for Kepler and ~ 1 kpc for CoRoT), which is small 
compared to the distance to the Galactic center. For our ex- 
ploratory analysis it is therefore safe to neglect the enhance- 
ment factor 77 variation within the volume probed by these 
stars. However, its dependence on the assumed strings ten- 
sion cannot be ignored as it spans a range from r\ ~ 10 3 ' 6 at 
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to rj » 10 5 - 10 5 ' 2 at G/x/c 2 < 111 1 ! :i . Plug- 



ging these numbers into |5| and integrating along the line of 
sight, one obtains 
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where the upper estimate corresponds to the survey parame- 
ters characteristic of Kepler (and high tension of strings with 
rather conservative consequences for r) and the lower is more 
relevant to CoRoT asteroseismology survey (with G/j,/c 2 close 
to the middle of the range we consider). 

The number of target stars is of order 10 5 for Kepler and 
10 — 100 for CoRoT so these estimates are not that small 
especially taking into account that fi s is poorly constrained 
at present. However, what is truly crucial here is the accuracy 
of observations given that the cosmic strings crossings are not 
recurring events 2 and any detection of strings would be one of 



Strings do oscillate but the time scale of these oscillations is of 
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those extraordinary claims that, by common wisdom, require 
extraordinary evidence. As a consequence, the light curve 
should be measured with high precision and well sampled to 
allow an unambiguous identification with an event described 
by pi. 



That is, we choose 



While the precision of photometry for both missions is very 
high, the sampling rate of Kepler is unlikely to be sufficient 
for our purpose. Kepler is loo king at Solar-type stars with 
R ~ 10 1 cm and according to (14 1 the duration of events is 
t ~ (10—20) s. Kepler integrates light for 30 min ~ 2- 10 3 s and 
therefore there is no chance to resolve the light curve of the 
string events, except for an unlikely case of string crossing 
the star almost in parallel to the string itself 3 . Thus, for 
Kepler targets x\ ~ and this mission has more or less no 
constraining power with regards to cosmic strings. However, 
when interpreting light curves of stars in search for planet 
transits one should keep an open eye to the possibility that 
an 'anti-transit', in which the brightness of the star suddenly 
increases according to ( 12 I could be caused by a cosmic string 



rather than a planetary transit. 

CoRo T case seems more promising 4 . In the asteroseismol- 
ogy mode of the satellite, the brightness of a small number of 
stars is sampled once a second (of which ~ 0.8 s is the integra- 
tion time and the rest is readout etc.), although the publicly 
available data are integrated to 32 s exposures, which is opti- 
mal for stellar seismology studies. The targets are relatively 
bright stars (my ~ 6 m — 9 m ) and despite the short exposure, 
their brightness is measured to a sufficiently high precision 
(from A/ « 0.005 for m v = 6 m to A/ « 0.02 for m v ~ 9 m ) 
at a reliable signal-to-noise of S/N = 10. This is roughly the 
range of parameters expected for cosmic string crossings as 
explained in the first two paragraphs of this subsection. We 
therefore proceed to see how a non-detection of such events 
in the CoRoT data constrains the parameters of loops. 

To do so, we need to specify what we would consider a 
reliable detection of the event - that is, impose specific 'cuts' 
on the parameters of events that would guarantee that it is 
noticed in the data and not confused with some other effect 
such as stellar variability. We choose the following criterion: 
(a) the duration of the event should be at least t m i n = 10 s 
so that the light curve ( 12 I is sampled at 10 points by the 



CoRoT satellite and (b) the light curve of the crossing event 
should be measured with an accuracy allowing to resolve it 
into at least K = 5 flux levels at a given signal-to-noise S/N. 



order cosmic string length divided by the speed of light, which is 
hundred and thousands of years. Moreover, their oscillations are 
not restricted to any particular mode and given the plentitude of 
available modes there unlikely to be any obvious periodicity in 
the recurrence of crossings; worse still, besides oscillations there 
is also center-of-mass motion of the string of unknown magnitude 
and direction. This makes recurrence of such events in a given 
source rather unlikely. 

3 Just how unlikely depends on which oscillations modes are ex- 
cited on the string. It appears that Pp = Px/iPx) i s a sensi- 
ble upper limit to the probability that transverse velocity of the 
string is less than j3j_. 

4 The CoRoT space mission was developed and is operated by the 
French space agency CNES, with participation of ESA's RSSD 
and Science Programmes, Austria, Belgium, Brazil, Germany, 
and Spain. 



1, t > t min and f max /K > (S/N)e F 
0, otherwise 
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where ef = cff/F is the observational relative uncertainty 
of the measured flux F, which is available from the CoRoT 
data. For a given source (R, D oa ) conditions in ( 18 1 limit the 
string tension from below at a value /x m i n that satisfies both 
the duration and flux increase thresholds. Assuming that 



the string is halfway to the source, one obtains, using ( 13 1 
and JIS 
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The efficiencies Xi are either 1 or depending on whether |19| 
is satisfied, thereby limiting our sensitivity at low tensions. 
Limitations at the other end of assumed tension range come 
from a gradual decrease of the enhancement factor r\ for 
Gn/c 2 ^ 10 -13 . We approximate the dependence found in [3] 
with a simple broken line 
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Table [I] presents (T,D os ,R,ef) for CoRoT target stars 
mined from NStED database [55] . Stars are included in the 
table and in the analysis regardless of their known variabil- 
ity. We assume that intrinsic stellar variations are unlikely 
to present themselves because string crossings with profile 
like (121 at these time scales; external variability is also ex- 



pected to affect the flux at much greater time scales, as the 
velocities of potential perturbers other than strings are not 
relativistic. 

Out 30 CoRoT targets, for which the data has been made 
public so far, only eight (top in the table) have reliable esti- 
mates of both distance and radius in the NStED database. For 
the majority (14 stars, middle of the table), only an estimate 
of distance is available while the radius has been estimated 
photometrically; this is of special concern in the case of giant 
stars (marked with * in the table) due to large uncertainty 
in estimating the bolometric correction for these stars. For 
eight stars (bottom of the table) both the distance and radius 
were estimated photometrically. We have checked numerically 
that giants and stars with photometrical distance estimates 
contribute little to the string density-constraining power of 
CoRoT data. 

Figurej3] shows the constraint on £} s as a function of fj, given 
by ( 10 111 assuming a non-detection of lensing signal in the 
light curves of stars listed in Table [I] One can see that the 
presently available CoRoT data does not significantly con- 
strain the density of cosmic string loops. This is mostly due 
to proximity of the source and strict cuts imposed on the lens- 
ing parameters, which severely limit the range of tension for 
which the lensing effect is assumed detectable. We, however, 
believe that these cuts are adequate given that cosmic strings 
have never been detected before and therefore strong evidence 
would be needed to lay such claims. 

The envelope curve of the constraints roughly follows oc 
^-3/2 ij ne w jt n occasional wiggles due to eventual inclusion 
of sources with larger and large /i m i n into analysis. This be- 
havior is due to the fact that the number density of strings, 
which controls the detection rate, is oc /i -1 at a fixed mass 
density, which is itself oc rj(fj,) oc pi~ 1//2 in the range of tensions 
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TABLE I: CoRoT target stars. The following parameters 
are given: the duration of the observing run T, distance to 
the target D, radius of the star R, flux measurement preci- 
sion ef- Where possible, we use D and R quoted in either 
NStED (preferably) [25] or SIMBAD database; where not, 
we compute these quantities from the quoted effective temper- 
ature, apparent and absolute magnitudes and color [55]. We 
estimated Ef as a square root of the average flux in photons 
because photon noise is dominant in CoRoT photometry of 
these bright stars. The last column shows the minimum ten- 
sion ( 19 1 probed by a given target star (assuming i m i n = 10 s, 
K = 5 and S/N = 10). 



Star id 


T (days) 


D (pc) R (R Q ) 


s F (%o) 


Gfi 


• /c 2 

nm/ ^ 


Golden: 
















HD 49933 


61 


OU 


l.o 


0.32 


2.0 


1U 


LI 


HD 51106 


61 


i on 
iyu 




0.68 


2.2 


i n~ 


L2 


HD 175272 


27 


or, 

O'J 


1 £ 
l.o 


0.70 


2.7 


1U 


L2 


HD 175543 


27 


i fin 


9 8 


0.59 


1.9 


i n~ 


12 


HD 175726 


27 


97 


i n 


0.48 


7.3 


1U 


LI 


HD 181420 


157 


1 

-±y 


1.0 


0.47 


2.8 


1U 


12 


HD 181440 


157 


1 «C\ 


a n 


0.29 


1.3 


1U 


L2 


HD 181906 


157 


71 


1 . 


0.75 


2.9 


i n~ 


L2 


HD 50846 


58 


1 1UU 


1 G 


1.12 


9.0 


1U 


L3 


HD 174884 


27 


you 


9 8 


0.91 


4.8 


1U 


L3 


HD 174966 


27 


1 9P1 


9 n 


0.81 


2.4 


1U 


L2 


HD 175869 


27 




b.y 


0.30 


1.1 


1U 


12 


HD 180642 


157 


480 


2.0 


1.03 


7.7 


10" 


L3 


HD 180973 


157 


100 


3.1 


0.49 


2.7 


10" 


12 


HD 181231 


157 


1100 


3.1 


1.24 


6.3 


10" 


L3 


HD 181555 


157 


110 


8.2 


0.72 


9.6 


10" 


12 


HD 182198 


157 


770 


4.8 


0.90 


1.0 


10" 


12 


HD 50170* 


61 


690 


160 


0.51 


2.1 


10" 


11 


HD 50747* 


61 


150 


23 


0.29 


7.9 


10" 


12 


HD 50890* 


55 


390 


410 


0.33 


6.2 


10" 


11 


HD 175679* 


27 


160 


12 


0.35 


4.8 


10" 


12 


HD 181907* 


157 


110 


14 


0.30 


6.8 


10" 


12 


HD 50773 


58 


330 


1.9 


1.82 


1.9 


10" 


12 


HD 50844 


58 


360 


2.7 


1.54 


2.1 


10" 


12 


HD 50405 


55 


330 


1.4 


1.70 


1.3 


10" 


12 


HD 292790 


55 


330 


3.0 


1.84 


3.0 


10" 


12 


HD 174936 


27 


220 


1.9 


1.23 


1.9 


10" 


12 


HD 175542 


27 


330 


2.1 


1.48 


1.7 


10" 


12 


HD 174987 


27 


330 


2.5 


1.50 


2.0 


10" 


12 


HD 181072 


157 


310 


2.0 


1.48 


1.7 


10" 


12 



two reasons. These small objects are crossed by cosmic strings 
in milliseconds and are often highly variable. This forces us to 
consider only the 'ultimate' light curve with a twofold increase 
in brightness as candidate lensing events and focus on the 
brightest of sources so that this signal is reliably detectable 
on top of the photon noise and intrinsic variability of the 
source at millisecond level. 

The object that best matches the above criteria is Scor- 
pion X-l, which is in fact the brightest persistent X-ray source 
in the sky not including the Sun. It is a low-mass X-ray binary 
at D « 2.8 kpc away, in which a neutron star accretes matter 
from a low-mass companion. The size of the X-ray emitting 
region of the source is believed to be just that of the neutron 
star and can be safely assumed to be 2R ^ 100 km, which 
corresponds to r ~ 10~ 15 in the sky. Observations of this 
source by the Proportional Counter Array (PCA) onboard 
Rossi X-ray Timing Explorer (RXTE) satellite with a total 
exposure time T = 3.32 • 10 s s were used by [27] to discover 
trans-neptunian objects in the Solar system by occultations 
they produce in the light curve of Sco X-l. The duration of 
such events is of order milliseconds, which happens to be very 
suitable for our purposes. 

The flux of Sco X-l observed by RXTE PCA is F w 
10 counts/s, which allows one to detect doubling of flux at a 
signal-to-noise level S/N = 10 on time scales t — (S/N) 2 / F = 
10 -3 s. The intrinsic variability is known to be Poissonian for 
Sco X-l and at these time scales is significantly below RXTE 
instrumental fluctuations (see 28 and figure 1 thereof in par- 
ticular). The time scale t = 10~ 3 s corresponds to the string 
crossing time for a source of 2R ~ 100 km and is less than 
that for all strings with tension Gfi/c 2 ^ G/i m i n /c 2 = 10~ 16 , 
which are also those for which doubling of the flux occurs. 
Therefore, we expect that a twofold increase in the flux of 
Sco X-l would be reliably detected had it happened during 
the observations. 

However, no such events were observed in the light curve 
(cf. figure 1 of 27). Application of (10 1 then yields the 



following constraints on the density of cosmic string based on 
the RXTE data for Sco X-l: 



= 5.6-10 



-2 ( Gfx/_ 
I 10 



" 1 Wri) 



(21) 



The constraints from RXTE observations of Sco X-l are plot- 
ted in figure [3] One can see that this data provides competi- 
tive upper limits on the density of loops of the lightest cosmic 
strings. 



considered. Therefore, the overall minimum of the graph is 
attained, ceteris paribus for those stars that have lowest /i m i n . 
The latter quantity is mostly affected by the angular size of 
the target star and this suggests to look for effect in the small- 
est sources possible. 



C. X-ray variability: Sco X-l 

Some of the smallest sources one can see in the sky are high- 
energy astrophysics sources such as accreting neutron stars or 
black holes. An object of R ~ (10 — 100) km placed at D ~ 
(1 — 10) kpc allows one to probe strings with tensions down 
to G/x/c 2 ~ 10~ 16 still in the small-source regime. However, 
using accreting neutron stars and the like is complicated for 



IV. SEARCH WITH PULSAR TIMING 

Usually pulsar timing (PT) is used to constrain the pres- 
ence of strings in rather indirect way: oscillations of string 
loops produce specific gravitational wave background and its 
influence can be sought in the anomalous residuals of millisec- 
ond pulsars. In this paper a different effect is considered. We 
look for anomalous residuals caused by the string crossing of 
the line of sight to the pulsar via the Kaiser-Stebbins effect. 

Such crossing instantaneously changes the apparent fre- 
quency of the pulsar by a small amount Av according to |2| . 
Observationally, this effect is the same as that of the gravita- 
tional wave burst with memory except that the latter applies 
to all observed pulsars whereas the former is specific to the 
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TABLE II: Pulsar timing array target pulsars including PPTA 
and results from Arecibo (*) and Arecibo+GBT(**). The 
following parameters are given: the rms of measured pulsar 
timing noise a, distance to the target D os |32) . total time span 
of observations T; we also list tension thresholds /i m i n probed 



Jl'M-2Mi (>ms - 0,152 jb) sost-'il 



by the pulsar as given by ( 25 1 



Pulsar name 


a (lis) 


D (kpc) T (yr) G/i m 


n/c 2 • 10 


J0437-4715 


0.2 


0.16 


4.3 


3.3 


J0613-0200 


1.1 


0.48 


5.5 


14 


J0711-6830 


1.6 


1.04 


4.4 


26 


J1022+1001 


2.2 


0.40 


5.5 


28 


J1024-0719 


1.3 


0.53 


5.5 


17 


J1045-4509 


3.0 


3.24 


5.2 


41 


J1600-3053 


1.0 


2.67 


5.5 


13 


J1603-7202 


1.9 


1.64 


5.5 


24 


J1643-1224 


1.7 


4.86 


5.4 


22 


J1713+0747 


0.5 


1.12 


5.5 


6.0 


J1730-2304 


1.9 


0.51 


4.6 


29 


J1732-5049 


3.5 


1.81 


5.5 


45 


J1744-1134 


0.8 


0.48 


5.5 


10 


J1824-2452 


1.7 


4.90 


3.1 


39 


J1857+0943* 


1.0 


0.91 


8 


12 


J1909-3744 


0.6 


1.14 


5.5 


8 


J1939+2134" 


~2.0 


17 


8.33 


8 


J2124-3358 


2.4 


0.25 


3.8 


45 


J2129-5721 


1.2 


2.55 


5.5 


15 


J2145-0750 


1.1 


0.50 


4.3 


18 



pulsar, whose line of sight was crossed by the string. Obser- 
vational signatures of gravitational wave bursts with memory 
were thoroughly studied in [29] and we will extensively utilize 
that analysis in the present work. 

The shift in frequency manifests itself in the pre-fit pulsar 
timing residuals as a broken line function that is equal to 
zero before the moment of crossing and grows linearly after 
that. Standard reduction procedures of pulsar timing data 
necessarily includes fitting of residuals for a priory unknown 
frequency and its first derivative, which absorb linear and 
quadratic trends correspondingly in post-fit residuals (see the 
middle panel of fig. j2j). 

We will assume that the string crossing is detected if the 
amplitude Ss of the residual due to the Kaiser-Stebbins effect 
exceeds twice the root mean square (rms) of the pulsar timing 
noise a: 



5s ~ 2cr, . 

Using eqs. (18-20) of we arrive at: 

3 Av r 



8s 



64 v 



-T. 



(22) 



(23) 



where T is the total time span of observations of the pulsar. 
This assumption corresponds to the following choice of x: 



1, s > 2a 
0, otherwise 



which is equivalent (cf. 23 1 to 



1, Gfi/c 2 ^ 64a/(37r 2 r) 
0, otherwise 



(24) 



(25) 



I J 



[ ■ , 




J17M-2M4 (rms = D.191 ps) past-'Tl 



■ 



-1000 

HJD-54B29.D 



FIG. 2: Pulsar timing residuals generated with the fake 
plugin of tempo2 package [30]. The first picture shows 10 
year span of pulsar timing with a — 150 ns, the second one 
present impact of string with G/i/c 2 = 10 -15 crossing at the 
fourth year of observations (timing noise omitted); the third 
one shows this impact on 'real' pulsar with a — 150 ns 
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FIG. 3: Combined constraints on the average density of cos- 
mic string loops based on CoRoT , RXTE and PPTA data at 
P=95% level. One can see that pulsar timing currently pro- 
vides the strongest constraints down to G/i/c 2 ~ 10~ 14 while 
lighter strings are constrained by the available data on Sco X- 
1; CoRoT data do not significantly constrain the population 
of loops at present. 



We are now ready to compute the constraints on the cosmic 
string loop density Q. s according to (jlO 11 1 based on the pul- 



according to |10 
sar timing data collected in the Parkes Pulsar Timing Array 
(PPTA) [31] project as well as from Arecibo [33] and Green 
Bank Telescope (GBT) [34] observations. The details of these 
data are listed in Table|H] The resulting constraints are shown 
in figure [3] which shows that current PT observations limit 
the density of light cosmic strings quite significantly. 



V. RESULTS AND DISCUSSION 

Figure [3] presents our final result including constrains from 
pulsar timing, X-ray data for Sco X-l and precision photom- 
etry from CoRoT . One can see that the most stringent lim- 
its come from pulsar timing except for the lightest strings 
where competitive constraints are provided by RXTE data. 
Presently available CoRoT data do not have much constrain- 
ing power in the range of tensions we consider. Existing obser- 
vations allow one to limit the average density of cosmic string 



loops down to Q. s ~ 10 at Gfi/c 



For larger ten- 



sions the limits become eventually weaker, proportionally to 
fj,- 1 for G/i/c 2 sC 10 -13 and /x" 1 ' 6 for Gfi/c 2 > 1(T 13 because 
of lesser enhancement of density of heavy strings. 

These results clearly demonstrate that the density enhance- 
ment in the Galaxy improves the limits on cosmic string abun- 
dance that can be set by test sensitive to the local population 
of loops as suggested by [3] . However, the same approach can 
also be applied to extragalactic sources, most notably quasars. 
The probability of lensing by string loops for these sources is 
lower due to absence of enhancement in the extragalactic case. 
On the other hand, the huge distance to the source can easily 
compensate for the relative deficit of string density. Moreover, 
these observations are insensitive to theoretical uncertainty in 
predicting rj. 



As an example we consider the 'Einstein cross' quasar 
QSO 2237+0305. The Optical Gravitational Lensing Experi- 
ment (OGLE) data for this source currently span more than 7 
years of high-quality photometry sampled every few days [35j 
and no obvious signature of string lensing are seen in these 
light curves 5 . This non-detection can be used to place some 
limits on the density of cosmic strings in the range of tensions 
3 • 10 -12 < Gfi/c 2 < 3 • 10~ 10 (corresponding to the strings 
with crossing time equal to the sampling interval and data 
span, respectively). 

Using | |11[ ) with enhancement factor T] = 1, D ~ 2 Gpc, 
T ~ 7.5 years yields 



4- 10" 



3- 10- 



(26) 



These limits are inferior to constraints from pulsar timing by 
about an order of magnitude. 

Future observational projects will noticeably improve the 
limits established in this paper. In J3B] it was proposed to use 
upcoming large-scale observational surveys of quasar variabil- 
ity to search for cosmic strings. In line with this paper, we 
might assume that a large number N ~ 10 3 of quasars is 
monitored for T ~ 10 years on a daily basis. Non-detection of 
string crossings in such survey would limit the average cosmic 
string loop density with Gfi/c 2 ^ 10 -13 ) at the level of 



n s = 2 • 10 



V 10- 13 



lGpc 
D 



10" 



10 yr 
T 



(27) 



which is orders of magnitudes better that constraints pre- 
sented in this paper - essentially as a by-product of AGN 
long-time variability study. 

Other avenues to better constraining f2 at lower tensions 
exist. With the eventual arrival of the Square Kilometer Ar- 
ray (SKA) we can expect further enhancement in the sensi- 
tivity due to the increase of number of pulsars and timing of 
more distant pulsars [37] . Pulsar timing array at SKA will 
consist of ~100 pulsars, that would be timed with precision 
better than 100 ns. Exact limits will depend on distances of 
pulsars, but we can forecast that accessible range and density 
estimates will improve by orders of magnitude. 
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The light curve for image B does show an approximately twofold 
increase in brightness starting near JD ~ 2453500 but present 
data is far from being sufficient to draw any firm conclusions. 
Unfortunately, no observations can resolve this source at present 
to the level needed to rule out the possibility of lensing by cosmic 
string. 
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We use the bolometric corrections from [38] . We compute 
distance 

£> = 10pcx io°-2( m v-Wv-A v ) 

In two cases, for stars HD 181231 and HD 182198, Ta- 
ble |T] quotes calculated values despite the presence of an 
estimate of D in NStED database. The reason for this is 
that the distance quoted by NStED has a very large un- 
certainty and significantly exceeds photometric estimate 
thereby making the latter conservative. 



